MOMENTS AND ^-COMMUTATORS OF NONCOMMUTATIVE 

RANDOM VECTORS 



KRZYSZTOF DROZDZEWICZ AND WOJCIECH MATYSIAK 



Abstract. A method for computing the mixed moments of (not necessarily 
commutative) random vectors from the first order moments, the g-commutators 
between the annihilation and creation operators, and the g-commutators be- 
tween the annihilation and preservation operators, is presented. The method 
is illustrated by a relevant characterization of g-Gaussian vectors. 



1. Introduction 

There has appeared recently a number of papers on the problem of identifying 
probability measures via the annihilation, preservation and creation operators as- 
sociated with them. Building on previous works of Accardi, Kuo and Stan (see e.g. 
[I] ;[!])) Stan and Whitaker jTT] showed that the moments of a probability mea- 
sure on M d can be recovered from two families of commutators, namely the com- 
mutators between the annihilation and creation operators, and the commutators 
between the annihilation and preservation operators, provided that the moments of 
the first order are known. The authors presented a concrete method for recovering 
the moments and applied it to identify one-dimensional normal distribution and 
one-dimensional Meixner distributions, postponing any multidimensional or non- 
commutative computations to some subsequent papers. In [10] . the foundations for 
defining the joint annihilation-preservation-creation decomposition of a finite fam- 
ily of not necessarily commutative random variables were laid. The analysis of the 
key example from [10], the Meixner random vector (X,Y) of class Ml, partially 
rests on a noncommutative extension (stated without a proof) of the method de- 
veloped in [IT] . (For some recent advances on Meixner classes of non-commutative 
generalized stochastic processes consult [4] and [5].) 

Our aim in this note is to provide more insight on the method when applied in 
dimension greater than 1. We present a direct multidimensional application of a 
slight generalization of the method from [11] . The generalization is twofold: we 
allow the analyzed random variables to be not necessarily commutative and we use 
g-commutators 

[A;B] q := AB - q ■ BA 

(A and B are some operators, q e [—1,1]). Next we show concretely how the 
generalized method works by solving a g-commutator problem connected with q- 
Gaussian random variables defined by Bozejko and Speicher in [6], where a re- 
markable q version of the Fock space was introduced. As special cases, we get 
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some characterizations of the multivariate normal distribution and of semicircular 
families. 

2. NONCOMMUTATIVE PROBABILITY BACKGROUND 

2.1. A joint annihilation-preservation-creation decomposition. We follow 
Stan [10] when presenting the material from this Section. Let (H, (•,•)) be a Hilbcrt 
space over R and let X^X^, ■ ■ ■ , Xd be symmetric, densely defined linear operators 
on H . Denote by A the unital algebra generated by X\, . . . , Xd (unital means that 
it is assumed that I £ A with I being the identity operator on H). We fix a unit 
vector cj) £ H belonging to the domain of any g from A (we assume its existence) 
and call it the vacuum vector. The elements of A will be random variables and the 
pair (A,((>) will be a probability space supported by H. For any g £ A, we will call 
E[g] := (g<j), (f>) the expectation of the random variable g. 

One can show that for any g £ A there exists a polynomial p in d noncommutative 
variables such that g = p(Xi, . . . ,Xd). After Stan [10] we introduce the following 
equivalence relation: 

Definition 2.1. Let (A,(j)) and (A',<f)') be two probability spaces (with expecta- 
tions E and E'), supported by two Hilbert spaces H and if'. We say that random 
vectors {X\, . . . , Xd) with Xi £ A and (X[, . . . , X' d ) with X[ £ A' are moment equal 
if for any polynomial p in d noncommutative variables 

E \p(X 1 , . . .,X d )] = E' \p(X[, . . .,X' d )} . 

For n £ NU {0}, define 

F n = {p(Xi, . . . , Xd)(f> : p - a polynomial of total degree less than or equal to n}. 

Clearly, F n is a finite dimensional (thus closed) subspace of H and Fq C Fi C . . . C 
H. Put G-i := {0}, Go := F and as G n take the orthogonal complement of F„_i 
into i 7 ^ for n £ N. Each G n is the homogenous chaos space of order n generated 
by Xx,.. . ,X d . Similarly, 

oo 

W := G n 

n=0 

is called the chaos space generated by X\, . . . ,Xd- It is easily seen that H is the 
closure of A<j> := {g<j) : g £ A} in H. 

Since X\, . . . ,Xd are symmetric operators, it is not difficult to prove (see [1]) 
that XiG n _L Gk for any i £ {1, ...,d}, n £ NU{0} and non-negative integer 
k {n — l,n, n+ 1}. It implies for any i £ the existence of three 

operators a~{i), a (i) and a + (i) such that 

(2.1) X l = a-(i)+a°(i)+a + (i) 

(the domain of the operators in (|2.1I) is A<f) . We will call them the annihilation, 
preservation and creation operators, respectively, what agrees with the fact that 
for i £ {1, . . . , d} and n G N U {0}, 

o-(i)(Gf n ) C o°(i)(G„) C G„ and o+(i)(G n ) C G n+1 . 

For the future reference, we now list some easily derived properties of the intro- 
duced operators. We shall be using the symmetry of a°(i) and the duality between 
o + (t)'s and a~(i)'s, particularly the identities: 

(2.2) (a o (*M0> = (&a°(*)^>, (a + (i)cf>, <f>) = (0, a -(i)^> , i = l,...,d. 
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We also will need the following consequences of the fact that <j> € Go : 

(2.3) a~(i)<p = a and a°(i)<j> = E[Xi]<p, i=l,...,d. 

Similarly, it is easily observed that if n G N and a : {1, . . . , n} — > {1, . . . , d}, then 

(2.4) a-(a(l)) . . . a~ (a(k))X a(k+1) . ..X^ = 
for k > [n/2\. 

Remark 2.2. If X%, . . . Xd are classical (commutative) random variables defined on 
the same probability space (O, J 7 , P) and having finite moments of any order, one 
can consider H = L 2 (r2,J r , P) and = 1 (the constant random variable, equal 
to 1), and regard X\,... ,Xd as multiplication operators on the space Al C H. 
Obviously, XiXj = XjXi for all 1 < i < j < d. 

2.2. Catalan sequences. If n G N, e = (e(l), . . . , e(2n)) G {1, -l} 2 ™ and k G 
{l,...,2n}, set r k = E-" fc e(j)- We say that £ is a Catalan sequence if r 2 „ > 0, 
T2n-i > 0, . . . , r 2 > 0, T\ — and let C(2n) be the set of such sequences. 

The Catalan sequences are important in the context of this note for the following 
reason. If Xi , . . . , Xm are random variables from a probability space (A, 4>) and 
e — (e(l), . . . , s(2n)) G {1, — l} 2n , then it is easily seen that 

^a £ «(l)...a e(2rl) (2n)0,0) =0 

unless e 6 C(2n) (here we adopt the convention that a +1 (fc) = a+(fc) and a _1 (fc) = 
a~(k) for any , thus a ±1 (k) are connected with X k via (|2.ip \ 

2.3. q-Fock spaces and g-Gaussian random variables. We briefly outline here 
the basics of Bozejko and Speicher [6] theory of q-Fock spaces. For a real Hilbert 
space Hq with complexification H, its g-Fock space F q (H) is defined as the closure 
of the algebraic sum C(f> © ® „ H® n in the scalar product 

n 

(fx (g>... ®f n ,gi ® .. -®g m ) q = 8 m . n l ? M • n (/j'SHj)) > 

7r£S„ 3=1 

where is the vacuum vector, 5„ are permutations of {1, . . . , n} and 

kl = #{(m) : * < 3,ir(i) > 

The range of q is [—1,1] - if q = — 1 or q = 1 one must first divide out by the 
null space, and then one obtains the usual antisymmetric (Fermion) and symmetric 
(Boson) Fock space. For q — the construction gives the full (Voiculescu) Fock 
space. 

Now, for f E H we define the creation operator a + (f), and the annihilation 
operator a~ (/) on C(f> (J) £f ®n m ^ e following way: 

« + (/)0 = /, 
a + (/)/i ®...®/„ = /(8>/i®...®/ n) 

and 

o~(/)^ = 0, 

n 

(/)/l ® • • • ® /n = 9,-1 ' ^> ' ^ ® • • • ® ® • • • ® /n- 

3=1 
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If q < 1 then a + (f) and a (f) extend to bounded operators on T q (H) and a (f) 
is (•, -) g -adjoint of a + (f). They satisfy the commutation relation [a~(f); a + (g)] q = 

if, 9)1. 

The self-adjoint operators x(f) = a + (f) + a (/) are called g-Gaussians. 

3. The main result 

3.1. A characterization. Having introduced the joint annihilation-preservation- 
creation decomposition, we are in a position to state the main result of the paper. 

Theorem 3.1. Let (A.,4>) and (B,ip) be probability spaces with the expectations 
E and E', supported by Hilbert spaces H and H' , respectively. Random vectors 
X = (X\, . . . , Xj) (with Xi s A) and Y = (Y%, . . . , Y^) (with Yi 6 B) are moment 
equal if and only if the following three conditions hold for all i,j £ {1, . . . , d}: 

(3- 1 ) [ a x( i )^ a x(j)] q = [oy(i);o^(i)] 5 ) 

(3-2) [<hc(*>>«x(j)] q = [o?(0;ay(j')] fl , 

(3.3) E[Xi] = E'[Yi\. 

(Here and later, the symbols a^{i), a^(i) and a^(i) will stand for the annihilation, 
preservation and creation operators associated with the i-th coordinate of a random 
vector U via (12.11) . ) 

Theorem 13.11 generalizes the main result of [Uj to g-commutators (note that 
there are no assumptions imposed on q - it can be an arbitrary real or complex 
number) and to not necessarily commutative random variables. 

3.2. Proof of the main result. 

3.2.1. Notation. Here we present some notation to be used in the remainder of the 
paper. 

While working with mixed moments of the form 

(3.4) E [X a{1) . . . X CT ( fe )] = (X a (i) . . . X ff (k)<t>, <f>) 

with fcsN and a mapping {1, . . . , k} into {1, . . . , d}, we will always put Zj = X a ^ 
for j = 1, . . . , k, so p. 41) will transform into (Z\ . . . Zk4>, <j>). 

Next, for the sake of brevity, we introduce the following notation: 

i:=Zi = X a(i) ,r := a^(o-(i)),i° := a n x (a(i)),i + := a\{a{i)); 

if A and B are any operators, then (A) will stand for (A<j), </>}, while the composition 
AB understood as (AB)x = A(Bx) will be denoted using comma, i.e. A, B := AB. 
For example, 

(l- 2- [3-; 4] q , 5,6) := (a x (a{l))a x (a(2)) [a x (a(3)); X a{4) ] g X a{5) X a{6) ^ 4>) ■ 

3.2.2. Proof of Theorem \S.l\ Essentially the proof will follow the same line of rea- 
soning as the proof of the main result of Stan and Whitaker [TT|. However, at some 
points we will need to take additional care of the issues arising from the lack of 
commutativity, and from the fact of using ^-commutation relations. 

The following Lemma generalizes " the product rule for commutators" |11( Lemma 
3.2] to q-commutators: 
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Lemma 3.2. For any operators Y, U x , . . . , U n 

n-l 

[Y; U x ... U n ] q = I" ■ Ui ...Ui [Y; U l+1 } q U l+2 . . . U n + (<f - q) ■ U x . . . U n Y. 

»=0 

Proof. We proceed by induction on n. It is easy to check that for any operators 
Y,Vy,Vi 

(3.5) [Y; V 1 V 2 ] q = [Y; V x ] q V 2 + q-V 1 [Y; V 2 ] q + (q 2 - q) ■ V X V 2 Y, 

so the result holds for n — 2. Assume it holds for some n > 2. Putting Vi = 
U\...U n and V% — U n+ \ in (|3.5j) we get 

[Y; C7i . . . U n+1 ] q = [Y;U 1 ... U n ] q C7 n+1 +g-C/i ...U n [Y; [/„ + i] g +(g 2 -g)-[/i • • • U n+1 Y. 
By the induction assumption 

n-l 

[Y; U x ... U n+1 ] q = <i ■ Ux ...Ui [Y; U i+1 ] q U l+2 . . . U n+1 + 

8=0 

(q n -q)-Ux... U n YU n+1 +q-Ux...U n [Y; U n+1 ] q + (q 2 - q) ■ Ux . . . U n+1 Y. 
Since Ux... U n YU n+1 = U x ...U n [Y; U n+1 ] q +q-U x ... U n+X Y, we obtain 

n 

[Y; U x ... U n+X ] q = <t-Ui --.Ui [Y; U i+l ] q U l+2 . . . [/„+i + (g" +1 - 9 )-C/i . . . U n+X Y, 

i=l 

and the proof is complete. □ 

Proof of Theorem \3.1\ Assume first that X and Y are moment equal. It is not diffi- 
cult to observe that the creation, preservation and annihilation operators associated 
with a random vector U are defined in terms of the orthogonal projections of the 
chaos space T-L generated by U x , . . . ,Ud on the relevant homogenous chaos spaces 
G n of order n (consult [T] for a detailed treatment of fully analogous commutative 
case). The projections in turn, from the Gram-Schmidt orthogonalizing procedure, 
are determined uniquely by the mixed moments of the vector U. Thus, from the 
fact that X and Y are moment equal it follows that their creation, preservation 
and annihilation operators are identical, so (|3.1[) and (|3.2j) hold. 

Assume now ||3TI]) . ([3T2| and ([3T3| . Fix k € N and cr(l), . . . , a(k) £ {l,...,d}. 
We will prove that 

E [X CT(1) . . . X a(k) ] = E' [Y a(x) . . . Y a(k) ] 

by showing that the expectation E[X CT /]\ . . . X a ^\ (and E'p^^) . . . Y^^]) can be 
expressed in terms of the expectations of monomials of lower degrees, and that the 
monomials depend only on the g-commutators between the annihilation and cre- 
ation operators, and the g-commutators between the annihilation and preservation 
operators. (Thus, by (|3.1I) and (I3.2|) the monomials will be identical when deter- 
mined for X using E and for Y using E'.) Repeating the process of degree lowering, 
we will eventually arrive at the expectations E[Xj] and E'[Yi] which are equal by 
(|3.3|) . This will end the inductive reasoning. 
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By ((23}, (HH and (l23)l . 

E . . . X CT ( fe) ] = (Zi . . . Z k (/), (/)) = 

<a^(o-(l))Z 2 . . . Z k <f>, 4) + (a° x {a(l))Z 2 . . . Z k cf>, 4>) + (a+{a(l))Z 2 . . . Z k , <j>) = 
(a x {a{l))Z 2 . . . Z k 4,4)) + (Z 2 . . . Z k ^,a\{a{l))4)) + (Z 2 . . . Z k ,a x {a{l))4^ = 

{a x (a(l))Z 2 . . . Z k <j>, 4>) + E [X a{1) ] ■ (Z 2 . . . Z k cp, 4>) . 

Since (Z 2 . . . Z k (f>, </>} is the expectation of a monomial of degree fc — 1 (in fact, the 
expectation of monomial X a ^ . . . X^^), we will focus on (a x (a{\))Z 2 . . . Z k (f>, 0). 
From (|2.3p and Lemma 13.21 it follows that 

(a x (a(l))Z 2 ...Z k <j>,<p) = (l-,2,.,.,k) = ( [l~; 2, . . . , k] q ) = 

k 

^ 2 • (2, . . . , » - 1, [I - ; i] , , » + 1, . . . , fe) + (q k - 2 - q) ■ (2, . . . , k, I") • 

Observe that the last term vanishes by another application of (|2.3p . Since the 
operator [l - ; i] q — [a x (cr(l)); X^^] q maps F k -i into F k -i, the expression 

(3.6) (2,...,*-l,[l-;*] gl * + l,...,fc) 

is the expectation of a polynomial of degree fc — 2 (for all i = 2, . . . , k). Thus we 
have expressed E[X ff (x) ■ • ■^■a(k)\ m terms of expectations of some monomials of 
lower degrees. 

What is left to show is that the monomials of lower degrees hidden in Q3.6|) depend 
only on the g-commutators [a^-(-); o^(')] and [a^-(-); o^(')] from assumptions 
fl3HJ) and (|3~2|) . To this end, we first expand via (|2~T|) : 

[1-;*],= [a x (a(l));X a(l) ] q = [a x (a(l));a x (a{i)) + a° x (a{i)) + a+(a{i))] q , 

so (a x (a(l))Z 2 . . . Z k cf>, (f)\ becomes the sum of three terms of the form 

k 

(3.7) <?~ 2 ■ (2, • • • , * - 1. [a x (a(l));a x (a(i))] q , i + 1, . . . , fc) 

i=2 

with e e {—,0,+}. Of course, the choices e = and e = + do not require any 
further arguments, so we focus on the case e = — . (The point is that we do not 
assume that [a x (-); a x (-)] q ~ [a Y (•); a Y (•)] q We have 

fc 

(3.8) ^ <T 2 • (2, . . . , i - 1, [a* (a (i))] q , i + 1, . . . , k) = 

i=2 

fc k 

^(7 l - 2 -(2,...,*~l,l-, l -, J + l,...,fc)-^ (? i - 1 .(2,...,z-l,*-,l-,* + l,...,fc). 

j=2 i=2 

Consider the first sum after the equality sign in (|3.8[) . Observe that there are 
fc — i operators Z that follow the rightmost annihilator in the product 2, . . . , i — 
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1,1 , i , i + 1, . . . , k. Lemma l3~2l and formula (|2.3|) imply that 

k 

Y / q i - 2 -(2,...,i-l,l-,r,i+l > ...,k) = 

fe-1 

q*- 2 ■ (2, . . • , i - 1, 1", [r ; i + 1, . . . , k] 9 ) = 

j=2 

fe-1 k 

E <T 2 E • (2, . . . ,i - 1, 1-, * + 1, . . . ,i - 1, [r ; j] g , j + 1, . . . , *) . 

t=2 j=i+l 

We use (|2.1[) again, and we again are interested only in 

(2, . . .,* - 1, + 1, . . . ,j - 1, [i~;j~] q , 3 + 1, • ■ • , fc) = 

(2,...,i- + ~ l,rj~J + l,...,fc) - 

(2, - l,l~,i + 1, - 1, + 1,..., k) . 

The number of the operators Z that follow the rightmost annihilator in the product 

2.. ...1-l,l-,i + l,...,j-l,i - ,j - ,j + l,...,fc (and2,...,i-l,l~,i + l,...,j- 

+ l,...,k) equals k — j, so it is strictly smaller than k — i for all j 
! + 1, . . , , k. Therefore repeating the above steps we will finally arrive at some 
products with the rightmost terms being some annihilators - the expectations of 
the products will vanish because of (|2.3|) . 

Clearly, the same reasoning can be applied to the second sum on the right-hand 
side of (pT8l) . 

Thus we have shown that the degree reduction procedure depends only on 
the g-commutators between the annihilation and creation operators, and the q- 
commutators between the annihilation and preservation operators, completing the 
induction. □ 

4. An example 

Let (^4, 4>) be a probability space with the expectation E supported by a Hilbert 
space H, and let q G [—1,1]. In order to present the algorithm described in the 
proof of Theorem 13.11 in action, we will solve the following q-commutator problem: 

Problem 4.1. Compute the moments of random vector X = (Xi, . . . , Xd) if Xi G 
A are such that 

[a-(i);a + (j)] q = c(i,j)-I, 

[o-(i);o°C7)], = 0, 
E[Xi] = 0, 

for i, j — 1, . . . , d, and the matrix (c(i, j))i,j=i ^ s positive definite. 
Some special cases of Problem l4.1l are discussed in Subsection 14.21 

4.1. A solution to Problem 14.11 
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4.1.1. Feynman diagrams. We now recall some basic facts on Feynman diagrams 
that are needed to solve Problem 14. 1( the notation we use is almost identical with 
the one from [7J Section 2], which should be consulted for more information. 

If S is a linearly ordered set with an even number of elements, then a complete 
Feynman diagram 7 on S is a partition of S into two-element sets. We shall denote 
the collection of all such diagrams by F c (5) and regard 7 as a set of ordered pairs 

(4-1) l = {(h,jl),--',(in,jn)} 

if S = {ii, i n , ji, ■ ■ ■ ,j n }, ii < ... < i„, ik < jk, jk / jl for fc ^ / (the j k are 
generally unordered) . We may represent a Feynman diagram with a simple graph, 
for example: 




1 2 3 4 5 6 7 8 



(this one corresponds to the Feynman diagram 7 = {(1,3), (2,7), (4,8), (5,6)} on 
{1,...,8}). 

A pair (k, I) £ 7 is a left crossing for G 7 if k < i < I < j and ci(i,j) will 
stand for the number of such left crossings for (i, j). The number 

c(7) = Cj ( i>J ') 

(*ij)&7 

will be referred to as the crossing number of 7 (restricted crossing number in 
Biane's [3] terminology). The crossing number can be determined by counting 
the intersection in the corresponding graph. For example, in the diagram above, 
c t (2, 7) = c(4, 8) = 1, q(1, 3) = cj(5, 6) = 0, and c( 7 ) = 2. 

For a given 2n-element set S and a Catalan sequence e G C(2n) we call a 
complete Feynman diagram on 5 of the form (|4.1[) compatible with e if e(ifc) = —1 
and e(jfc) = 1 for all k. If A C C(2n) then we denote by F C (S,A) the collection of 
all complete Feynman diagrams on S compatible with some eei. 

4.1.2. Some low-order moment computations. To get a feeling of how the algorithm 
described in the proof of Theorem 13 . 1 1 works under the assumptions of Problcm l4.11 
we will compute a few low-order mixed moments of random variables X\, . . . , Xd- 
The first moments E[Xj] are assumed to be zero. Let us compute the "covari- 
ances" : 

E,[XiX s ] = (ZiZaM) = (1> 2 ) = (i^ 2 ) = (t 1 ^ 2 ],) = 

<[i-; 2 1 9 ) + ([i-; 2 °] 9 } + ([i-; 2+ ] 9 }. 

The first term vanishes because (a~(i)a~(j) — q ■ a~(j)a~(i)) <f> — 0. From the 
assumptions of Problem 1431 it follows that ^[l-;2°]^ = and /[1~;2+]^ = 

([a-(i);a+(j)] q ) = c(i,j) (70,0). Hence E^X,} = c(i,j). 
The third moments are zero. Indeed, Lemma 13.21 implies that 

(ZiZ 2 Zs<t>,<l>) = (1,2,3) = <l-,2,3> = ([l-;2,3] g ) = ([l _ ;2] fl ,^+q-(2, [l-;3] g ), 
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and from the assumptions of Problem 14.11 and formulas (|2.4[) and (|2.3|) . it follows 
that both 



and 



( [1-; 2] q , 3) = c(a(l), a(2)) • E[X a(3) ] + 2~,3) - g • <2~ 1",3) , 
(2, [1~;3] 9 ) = c(a(l), a(3)) • E[X ff(2) ] + (2, 1~, 3") - q • (2, 3", 1~) 



vanish. 

The last example we present is a computation of the fourth mixed moments. We 
have 



(Zi . . . Z 4 0, 4>) = (1, 2, 3, 4) = (I" 2, 3, 4) = ( [I"; 2, 3, 4] q ) = 

' [I- ; 2] 3, 4) + q ■ (2, [1- ; 3] 4) + q 2 ■ (2, 3, [l"; 4] , 



Now the key point is the following. It turns out that a smart way of handling 
the three elements of the sum above is to expand only the first one via (|2.1[) (if 
one expands at each stage all three elements of relevant sums, the total number 
of elements under consideration grows fast and it is hard to find any pattern or 
regularity). Let us expand then: 



[l~;2],,3,4) = ([l-;2-] g ,3,4)+([l-;2 ] 9 ,3,4) + ([l-;2+]^3,4) = 

2" 3, 4) - q ■ (2-, 1-3, 4) + c(<r(l), a(2)) ■ (3, 4) . 

The term (3, 4) = (3 _ , 4) is the expectation of a monomial of a lower degree (in fact 
it equals c(er(3), ct(4)); we leave (1 _ , 2 _ , 3, 4) temporarily unexpanded and focus on 



-«-(2- l-,3,4) = -q-{2-, [l-;3,4] fl ) = -«-(2~, [l"; 3] q , ^-q 2 -^ , 3, [l";4] ( 
Thus we arrive at 

(1, 2, 3, 4) = c(ct(1), <t(2)) • c(<t(3), cr(4)) + 2", 3, 4) + 
?-(( 2 « t 1 ^ 3 ],' 4 ) - ( 2 ~' [ 1_ ; 3 ] g > 4 ))+? 2 -((2,3, [1";4] 9 ) - (2-, 3, [l-;4] g ; 
By ([23]) and ([231) applied to Z 2 in <^2, [1"; 3] g , 4^> and (2, 3, [1"; A] q 

(2, [l-;3] 9 ,4) = (2-, [l-;3] g ,4) , (2,3, [l";4] ? ) = (2-, 3, [l";4] 9 ) , 
so we are left with a computation of (1 — , 2~, 3, 4): 

(4.2) (1-2-3,4) = (l-, [2";3,4] 9 ) = (l~ [2"; Z] q , 4) + q ■ 3, [2-4], 
We again expand only Z3 in [2 _ ; 3] and obtain 

(1" 2~ 3,4)=c(ff(2)X3))-(l-,4) + (l-,2- 3~ 4)- 
g•(l-,3- 2~ 4)+g.(l- 3, [2";4] 9 ] 
Formula (H3) implies that (1~, 2~,3-,4) = (1",3-, 2", 4) = and 

'l-, 3, [2-; 4] \ = (l-, 3, [2";4+] \ = c(a(2), (7(4)) • (l~, 3) . 
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Concluding, we arrived at 

(1, 2, 3, 4) = c(a(l), a(2)) ■ c(a(3),a(A)) + c(<r(2), a(3)) • c(a(l), a(4))+ 

g- ^(2)^(4)). c(<r(l),<7(3)). 

4.1.3. XTie solution. Now, for the rest of Subsection l4.ll fix k 6 N and er(l), . . . , cr(fc) G 
{1, . . . , d}. We will be considering K = {ki, . . . , fcjv} (iV is a positive integer), sat- 
isfying 

(4.3) ki < . . . < k N and fej € {1, . . . , k}, i = l,...,N. 

If N is an even number (say N = 2n for some n G N), then for any K satisfying 
(|4.3j) and a Feynman diagram 7 = {(«i, ji), . . . , (i n , jn)} on if, we define 

n 

= n c ( cr ( fc ^)' cr ( fc i m )) ■ 

m— 1 

In order to solve Problem 14.11 we will compute the mixed moments of the form 

E [X a{1) . ..X a(k) ] = (X a{1) . . .X a(k) (f>, 4>) = {Zx... Z k cj),<i)) = (1, ...,k) . 

Our task will be to show that they match the ones of (/-Gaussian variables that are 
given by 

Theorem 4.2 (q-Wick theorem). For any q-Gaussian random variables Yi with 
"covariances" E[YjYy] = c{i,j), we have 

(4.4) E[Y a{1) ...Y a(2n) }= 9 C(7 M7) 

7£F c ({l,...,2n}) 

where 

n 

!/( 7 ) = ^({(itji), (i n ,j n )}) = Y[ c(cr(i m ),a(j m )) , 

m— 1 

and E[r tr(1 ) . . . Y^ 2n +i)] = 0- 

Theorem 14.21 was found by Bozejko and Speicher |6] (see also [7] for an alterna- 
tive treatment). It generalizes the famous Wick formula (see Subsection 14. 2\i for 
computing higher-order moments of the multivariate normal distribution in terms 
of its covariance matrix, which had apparently been discovered first by Isserlis [8] 
and rediscovered independently by Wick |12j . 

We begin with a lemma that captures the cancellation mechanism that has al- 
ready appeared in the computation of the fourth moments in Subsection 14. 1 .21 

Lemma 4.3. Fix N 6 N and K = {k%, . . . ,kjsr} satisfying (|4.3|) and define for 
1 < I < m < N 

R(l;m;U) := (fcf , . . . , k m , U) - (fcf , . . . , fcf , k m , U) , 

U is any operator. 

Under the assumptions of Problem \4-l\ the following statements hold: 
1 

(4.5) R(l; m;U) = J2 J~ j ' c a(k m )) ■ (fcf, . . . , kj , . . . , kf , u) 

i=i 
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(here and later the hat sign stands for the deletion of the relevant operator from the 
product), 

(4.6) (ki,...,k~,k m+ i,...,k N ) = 

m ^ 

^q m ~ J ■ c(a(kj),a(k m+ i)) ■ ^fcf, . . . , k~ , . . . , k~, k rn+2 , . . . , k N ^ + 
i=i 

(k-^ , . . . , k m+1 , k m +2, ■ ■ ■ j ftjv) ■ 
Proof. We start with a justification of the first statement. By Lemma [3721 
R(l;m;U) = (fcf , . . . , fcf ,k m ,U) - (fcf , . . . , hj~, k~, U) = 

(k x , . . . , fc^n [fc; ; &m]g 5 C/^ — (^i i • • • > [&J i fern] g ' ^) 

+ 9 ' {(^1 ' • ' ' ' ^-D ^mj [&J ' ^1g) _ (^1 ' ' ' - ' Vl'^mi [kj i ^1g) } 

From (|2.1|) it follows that the first difference equals 

c(a(ki),a(k m )) ■ (fcf , . . . , k^_ v U) ; 

the difference in the curly bracket is R{1 — l;m; [^~;C/] g ), which in fact equals 
R(l — 1; to; kf, U). Hence 

(4.7) R(l; m; C/) = c(a(h), a(k m )) ■ (fcf. • • • , V-i> U) + q ■ R(l - 1; m; fcf , IT). 
On the other hand 

i?(l; to; CO = <*r, km. U) - (K,k-, U) = c(a(h),a(k m )) ■ (U) + 

q-{(km,[KiU] q )-(k-,[kp,U] q )}. 

By (12.1[) . (|2.2[) and (|2.3[) . the difference in the curly brackets vanishes, so i?(l; to; U) = 
c( cr(fci ), <x(fc m ) • (C/). Combining (|4.7I) with this, we get 

(-1 

R(l;m;U) = ^2q 3 ■ c(a(k^ 3 ),a(k m )) ■ /fef, . . . , kf_ p . . . , kf, . 
3=0 

A simple reordering gives (|4.5|) . 

Formula (|4.6|) readily follows from the definition of R(m; to + 1; k n +2, • ■ • , fcjv) 
and (|4~5j) . 

□ 

Corollary 4.4. If n E N and K = {fc l7 . . . , fc 2 „_i} satisfies (|4.3p . then for all 
to 6 {1, . . . ,2n- 1} 

^fc-^ , . . . , fc m , fe m _j_i, . . . , fc2n — l) = 0. 

Proof. It immediately follows from (|4.6[) and the assumption E[X,-] = 0, j = 
1, . . . , d, combined with (|2.4[) . □ 

Corollary 4.5. For any m, n € N, m < n and K = {fc 1; . . . , fc2n} satisfying ()4.3[) . 

(4.8) (fcf ,...,fc m ,fc m+1 ,...,fc 2 „) = 

n+l i— 1 ^ 
i=m+lj = l 
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Proof. If m < n then iterating (|4.6[) one gets 

(ki , . . . , k m , fc m +l, . . . , fen) 

■ l-i ^ 

^2 ^Zl 113 • C ( Cr ( fc j)' fJ ( fc 0) • (fcr'---' fc 7:---: fc i-l' fc i+l:---' fc 2™) + 

(=m+lj=l 

^fe , . . . , k n , fe+1, • ■ • , fen) ■ 

On the other hand, from (|4.6p and (|2.4p it follows that 

(4.9) (fcf , ...,fc~,fe + i,...,fe„) = 

n ^ 
^g n ~ j ■ c(a(kj),a(k n+1 )) ■ (k±, ...,k~,.. .,k~,k n+2 , ■ ■ ■ ,k 2n ) ■ 
j=i 

□ 

As it was seen in the examples of low-order moment computations from sub- 
section 14. 1 .21 the algorithm from Theorem 13.11 applied to the expressions of the 
form 

(k l , • • • , k m , k m +\, . . . , fe n ) 
expands them, i.e. it transforms the operators Zk m+1 , . . . , Zk 2n either into the 
relevant annihilators, or creators (the expressions in which any of Z's is transformed 
into the preservation operator, vanish due to the assumption [a - (•); = 0)- 
Proposition ^. 6l stated below links the expressions with the Feynman diagrams: the 
annihilators in the expanded versions arc represented by the first elements of the 
ordered pairs from (|4.ip (the ascending lines in the corresponding graph) , while the 
creators are represented by the second elements (the descending lines in the graph). 

Proposition 4.6. For any m,n 6 N, m < n and K — {fe,...,fe„} satisfying 

631), 

(4.10) (fe-,...,fc-,fc m+1 ,...,fe n ) = J2 9 C(7 M7), 

7 6F c (_fs:,A m , 2 „) 

where A w] = {e = (e(l), . . . , e(j)) e C(j) : e(l) = . . . = e(i) = -1}. 

Proof. We proceed by induction on n. Clearly, (14.10[) holds for n — 1. Assume then 
it holds for some n-leN and all m < n — 1. 

Now, let m G {1, . . . , n}. By the induction assumption applied to (|4.8[) . 

^fe , . . . , k m , fc m _|_i , . . . , fen) = 
n+1 ;-i 

]T ^-z'-^-cKfe),^))- £ ^-"W- 

;=m+lj = l teFa(K\{k ] ,k l },A l - 2 ,2 n -2) 

For / 6 {m + 1, . . . , n + 1} and j G {1, . . . , I — 1}, consider any Feynman diagram 
7 from F c (K\{kj,ki}, Ai_i t 2 n --2) . If one extends 7 to the unique Feynman diagram 
7' G F c (if, Aj-i^n) by putting fe and A;/ into their natural places and linking them 
together, then = 1/(7) • c(a(kj),a(ki j) and 0(7') = 0(7) + Z — 1 — j, as there 

are I — 1 — j elements of K between fe and ki and none of the elements that are 
between fe and ki is linked with any other from these elements. (That is, there 
do not exist k a < fe, both from {fe+i, . . . , fe-i}, such that (fe, fe) G 7' - if they 
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existed, one would not have 7' E F c (K,Ai-i t 2n) because of e(/3) = 1. So ascending 
lines that correspond to all the I — 1 — j elements cross the lines joining kj and ki.) 
Hence 

n+l l-l 

(K,...,k-,k m+1 ,...,k 2 n) = ^2^2l c(j,) ' V W) 

with 7"s running through the set of all complete Feynman diagrams on K with kj 
and ki linked together, i.e. such that (kj,ki) G 7. 

Since all the 7"s from the formula above are in fact all the Feynman diagrams 
on K compatible with A mt 2n, we get (14. 10)) with m e {1, . . . , n}. 

□ 

Solution of Problem \4-l\ If k is an odd number then 

E[X a(1) ...X a ( k) ] = (l,...,fc) = (l-,2,...,fc) = 
by Corollary 14.41 If k = 2n for some n € N then 

E[X a(1) ...X a{k) ] =<l,...,2n) = (l- 2,...,2n)= ^ g c(7) • 

7eF c ({l,...,2n},Ai >2 „) 



the second equality follows from (|2.2p and the assumptions of the Theorem, and 
the third one is a consequence of Proposition ^. 61 Since 



F C ({1, . . . , 2n}, A h2n ) = F C ({1, . . . , 2n}, C(2n)) = F C ({1, . . . , 2n}), 
we finally arrive at the following formula: 

r -, 1 if A; is an odd number, 

lE 7 6F ({l,...,3n}) ? ( J ' "(7) lf fc = 

□ 

Thus, by Theorem 14. 2[ we get the following 

Corollary 4.7. Random variables X%, . . . , /rom Problem \4.1\ are moment equal 
to q-Gaussian random variables. In fact, by Theorem ] 3. 11 the assumptions of Prob- 
lem\4-l\ characterize random variables that are moment equal to q-Gaussians. 



4.2. Concluding remarks. Problem 14.11 covers some important distributions as 
its special cases. 

If q = 1 then we have the usual commutator [A, B] = AB — BA and g-Gaussian 
random variables may be identified with jointly Gaussian random variables. Un- 
surprisingly, from the solution of Problem 14.11 we get the following 

Corollary 4.8. If Xi € A are such that for all i,j £ {1, . . . , d} 

[a~(i);a + (j)] = c(i,j)-I, 
[a-(i);a°(j)] = 0, 
E[Xi] = 0, 

then 

E[X <7(1) ...X a{2n) ]= K7), 

7 eF c ({l,...,2n}) 

and E [X a{ i) ■ ■ ■ X^n+l)] = 0. 
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This is the celebrated Wick's formula, which is a special case of Theorem 14.21 
Since Gaussian distributions are determined by moments, the assumptions of Corol- 
lary characterize zero mean <i-dimensional normal vectors (in particular, when 
d = 1 we get the solution of Example 4.1 from 

Another interesting case is q — with [A, B]q = AB, and with such convention 
in the g-Wick Theorem that q c — for c ^ 0, and q° — 1. A Feynman diagram 7 
is noncrossing if 0(7) = 0. If NC(S') denotes the noncrossing diagrams on a set S, 
then 0-version of the g-Wick formula reads: if k = 2n then 

(4.11) E[x a(1) ...x (T(k) ]= Yl "M; 

7£NC({l,...,2n}) 

if k is odd then the expectation vanishes. Since (|4.11l) defines a semicircular family 
of covariance (c(i, j))ij=x,...,d (see [9| Definition 8.15]), the solution of Problem |4. II 
gives the following result: 

Corollary 4.9. If Xi € A are such that for all i,j £ {1, . . . , d} 

[a-{i);a + (j)] = c(i,j)-I, 
[a-{i);a°(j)} = 0, 
E[Xi] = 0, 

then they are moment equal to a semicircular family of covariance (c(i, j))i,j=i,...,d- 

By Theorem 13. 1[ both Corollaries are in fact characterizations. When q = — 1 
one gets a similar characterization of the Fermi field (which leads to ±l-valued 
symmetric random variables in classical probability). 
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